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We report site-selective oxygen isotope ( 16 0/ 18 0) effect (OIE) measurements on the in-plane 
penetration depth A a j, in underdoped Yo.6Pro.4Ba2Cu307_i, using the muon-spin rotation (/iSR) 
technique. A pronounced OIE on the transition temperature T c as well as on A~ 2 (0) was observed, 
which mainly arises from the oxygen sites within the superconducting Cu02 planes (100 % within er- 
ror bar). The values of the corresponding relative isotope shifts were found to be T c /T c = —3.7(4) % 
and AA~ i) 2 (0)/A~ b 2 (0) = —6.2(1.0)%. Our results imply that in this compound the phonon modes 
involving the movement of planar oxygen are dominantly coupled to the electrons. 



PACS numbers: 76.75.-K, 74.72.-h, 82.20.Tr, 74.25. Kc 

More than 15 years after the discovery of the cuprate 
high T c superconductors (HTSC), their pairing mecha- 
nism is still not understood. For a long time, a role of 
phonons in this mechanism was dismissed, despite evi- 
dence of their relevance, for example by isotope exper- 
iments (see, e.g., Refs. O ne 01 the reasons 

for this was the prevailing belief that phonons simply 
could not account for such high transition temperatures 
(T c ). The recent discovery of MgB2^ where phonons are 
clearly the pairing glue leading to an unexpectedly high 
T c Kt 40 K, clearly falsified this belief. Many other widely 
cited reasons as to why phonons ought to play no role in 
HTSC also do not stand to scrutiny^*! Recent work by 
inelastic neutron scattering (INS) 10,11 and angle resolved 
photoemission spectroscopy (ARPES) 12 ' 13 gave good in- 
dications of a strong electron-phonon coupling in various 
HTSC, invigorating the debate about the relevance of 
phonons i 1 ^* 1 ^ 

For conventional superconductors (including doped 
fullerensi^ and MgBzii), key experimental evidence for a 
phonon mediated pairing was provided by measurements 
of the isotope effect (IE) on T c . In complex materials, 
such as the HTSC, T c may depend on many parameters in 
a non-trivial way. Important complimentary information 
can be obtained from the IE on the magnetic penetra- 
tion depth A. Previous studie a 3 i 4 i 5 i 6 showed substantial 
oxygen isotope ( 16 0/ 18 0) effect (OIE) on the in-plane 
penetration depth A Q b, which indicates an unusual (e.g., 
non-adiabatic) coupling of the electrons to phonon modes 
involving the movement of the isotope substituted atoms. 
It is important to identify the relevant phonon modes re- 
sponsible for this effect. Keen insight in this respect can 
be obtained from studying the site-selective oxygen iso- 
tope effect (SOIE), as we will detail below. While two 
reliable reports on the SOIE on T c exist^J^ the SOIE on 
X a b was previously deduced only in a very indirect way 
from magnetization measurements in the vicinity of T c i& 



In this communication, we report on the SOIE on T c 
and on A a t in underdoped Yo.ePrcuE^CusOy-a, mea- 
sured directly and accurately by using the transverse 
field muon spin rotation (TF-/iSR) technique^ The re- 
sult clearly indicates a strong coupling of the electronic 
subsystem to phonon modes involving movements of the 
oxygen atoms in the Cu02 plane, while suggesting that 
modes involving apical and chain oxygen are less strongly 
coupled to the electrons. This agrees with the con- 
clusions observed from the combination of INSU2*ii and 
ARPESiSii 3 . results, which were obtained in a completely 
different way and in other HTSC. 

Polycrystalline samples of Y .6Pro.4Ba2Cu307_5 were 
prepared by standard solid state reaction^ In order 
to prepare oxygen site-selective samples a two-step ex- 
change process was appliediS* 2 ^ In the first step [500°C, 
35 h at 1.2 bar in 16 02 ( 18 02) gas] two pairs of completely 
oxygen substituted samples ( 16 O pac and 18 O pac ) were 
prepared [here indexes p, a and c mean planar (within 
Cu02 planes), apical and chain oxygen, respectively]. In 
the second step (330°C, 150 h, 1.2 bar) two site-selective 
oxygen samples ( 16 O p 18 O ac and 18 O p 16 O ac ) were pre- 
pared via annealing one 16 O pac sample in a 18 02 atmo- 
sphere and one 18 O pac sample in 16 02 gas. The other 
two samples (one 16 O pac and one 18 O pac ) were simulta- 
neously annealed in the same atmosphere as before. For 
the totally exchanged sample ( 18 O pac ) the 18 O content, 
estimated from the change of the sample weight after an- 
nealing, was found to be 90(2) %. The oxygen-tracer 
diffusion calculation^ for the material used in this work 
(average measured grain size is 16 fim) shows that after 
annealing at 330°C, during 150 h more than 90 % of the 
apex+chain sites and only about 10-15 % of the plane 
sites should be exchanged. Keeping in mind that there 
are 4 planar and 3 apical and chain oxygen sites in the 
unit cell of Yo.6Pro.4Ba2Cu307, for the above mentioned 
annealing procedure, one would expect about 45 % of an 
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average oxygen substitution. The experimental estimate 
(from the change of the sample weight) gives 41(2) % for 
both site-selective samples, in good agreement with the 
theoretical calculation. 
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FIG. 1: Room-temperature Raman spectra of the 
Yo.6Pro.4Ba2Cu307_5 samples. For xx polarization the line 
corresponds to the out-of-phase motion of the planar oxygen 
atoms ( 16 0, 325 cm -1 ; ls O, 294 cm _1 ). For zz polarization 
the two lines correspond to the in-phase motion of the CuC>2 
plane oxygen ( 16 0, 436 cm -1 ; ls O, 415 cm -1 ) and to the 
bond-stretching mode of the apical oxygen ( le O, 512 cm" 1 ; 
ls O, 482 cm" 1 ). 

The site-selectivity of the oxygen exchange was then 
checked by Raman spectroscopy. The 488 nm and 
514.5 nm lines of an Ar laser were used for excitation. 
Figure^shows Raman spectra of all samples with zz and 
xx polarizations. In the zz polarization for the 16 O pac 
sample we observe two lines, one corresponding to the in- 
phase motion along the c-axis of the planar oxygen atoms 
(436 cm -1 ) and one corresponding to the bond-stretching 
mode of apical oxygen, which in this compound is located 
at 512 cm -1 i2i In the xx polarization a line at 325 cm -1 
is seen, corresponding to the out-of-phase vibrations of 
the planar oxygen atomsiSl In the 18 O pac sample, the 
Raman lines are all shifted to lower frequencies with 
amounts agreeing well with the ones obtained in Ref. |2, 
indicating a nearly complete exchange of 16 O with 18 O. 
In the site-selective sample 16 O p 18 O ac , only the position 
of the apical oxygen line is shifted to the lower frequency, 
whereas the lines corresponding to the plane oxygen stay 
the same [apart from a small shift of one Raman line 
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FIG. 2: Section near T c of the normalized (to the value at 
5K) magnetization curves (1 mT, FC) of the site-selective 
Yo.6Pro.4Ba2Cu307_5 samples. 



(433 cm -1 instead of 436 cm -1 ), probably due to a small 
unintentional partial substitution by 18 O]. Note that the 
apical line (492 cm -1 ) is also slightly shifted from the 
expected 482 cm -1 , indicating that oxygen exchange for 
the apical and chain oxygen is slightly incomplete. In the 
18 O p 16 O ac sample only two planar lines are shifted, while 
the apical line stays the same. Therefore, the Raman re- 
sults confirm the almost complete site-selective oxygen 
substitution in all samples. 

To determine the site-selective oxygen isotope shift on 
T c , low field [lmT, field cooled (FC)] SQUID magnetiza- 
tion (M) measurements were performed. Figure [5] shows 
the M(T) dependencies in the vicinity of T c for all four 
samples. It can be seen that the 16 O pac sample has a 
higher T c than the 18 O pac one with nearly the same tran- 
sition width. The isotope shift AT C = -1.66(17) K for 
completely substituted samples ( 16 O pac and 18 O pac ) is 
in agreement with previous results**^ (T c was extracted 
from the intercept of the linear extrapolation as shown in 
Fig. [21). The data points of the site-selective 16 O p 18 O ac 
sample coincide within error with those of the 16 Dac 
sample, whereas the 18 O p 16 O ac data points coincide with 
the 18 O pac ones, indicating that the OIE on T c is related 
to the oxygen within the CuC>2 planes. The results of 
the SOIE on T c are summarized in Table U The OIE ex- 
ponent on T c was determined as qo — — d In T c /d In Mo 
(Mo is the mass of the oxygen isotope) . 

In order to determine the SOIE on A a fc we performed 
TF-//SR measurements at the Paul Scherrer Institute 
(PSI), Switzerland, using the 7rM3 beam line. The sam- 
ples were cooled from far above T c in a magnetic field of 
200 mT. In a powder sample the magnetic penetration 
depth A can be extracted from the muon-spin depolar- 
ization rate a(T) oc 1/A 2 (T), which probes the second 
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into evidence these results, the curve fitting the tempera- 
ture dependence of a sc {T) for 16 O pac (shown below) was 
subtracted from the experimental data [Fig.^b)]. It can 
be seen that the experimental points for the two pairs of 
samples mentioned above coincide within the error bars 
(nearly within the error bars in the case of 16 O pac and 
16 O p 18 O ac ). The data in Fig. 0| were fitted to the power 
law a sc (T)/a sc (0) = 1 - [T/T c ) n with cr sc (0) as a free 
parameter^ The value of T c was taken from the mag- 
netization measurements (see Table 0, n was fitted as a 
common parameter for all four data sets and was found to 
be n = 1.87(4), which is typical for underdoped YBCO. 19 
Note, that the exponent n is isotope independents^ and 
that fits with the same n for all data sets allow to in- 
crease the accuracy in the determination of er sc (0). The 
values of cr sc (0) obtained from the fits are listed in Ta- 
ble [I] and are in agreement with previous results^* 2 ^ 
From the values of a sc (Q) the relative isotope shift 
of the in-plane penetration depth AA ah 2 (0)/A afc 2 (0) = 



,)(0) and the OIE 



(X)(0)-a sc ( lb O pac ){0)}/a sc ( lb O pacJ 

din. AT. (0) /din Mo were determined 
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exponent flo = 
(see Table HJ) 

and 18 O pac . Our results (summarized in Table are in 
good agreement with previous estimates of SOIE from 
magnetization experiments^ and our previous /iSR non 
site-selective OIE results^ 



FIG. 3: (a) Temperature dependence of the depolariza- 
tion rate a sc in site-selective Yo.6Pro.4Ba2Cu307_a samples 
(200 mT, FC). Data points below 10 K are not shown (see 
text for an explanation). The solid lines correspond to fits 
to the power law a sc (T)/a sc (0) = 1 - (T/T c ) n for the 16 O pac 
and 18 O pac samples, (b) Same as in (a), but after subtracting 
the fitted curve for the ie O P ac sample. 



moment (AS 2 ) 1 / 2 of the probability distribution of the 
local magnetic field function p{B) in the mixed stated 
For highly anisotropic layered superconductors (like the 
cuprate superconductors) A is mainly determinated by 
the in-plane penetration depth A a h : cr(T) oc l/A 2 b (T). 19 
The depolarization rate a was extracted from the /iSR 
time spectra using a Gaussian relaxation function R(t) — 
exp[— a 2 t 2 /2]. As discussed in Ref. 0, the measured to- 
tal depolarization rate a is determined by three terms: 
a superconducting (cr sc ), an antiferromagnetic (cr a / TO ), 
and a small nuclear magnetic dipole (a nm ) contribution. 
The antiferromagnetic contribution o~ a f m is present only 
at low temperatures. Therefore data points below 10 K 
were excluded in the analysis. The dipole contribution 
o~nm is practically temperature independent below 100 K, 
and o~ sc was determined by subtracting o~ nm measured 
above T c from a according to: a sc = \J a 1 — a 2 lm . Fig- 
ure |31 shows the temperature dependence of a sc for the 
Yo.6P r o.4Ba2Cu307_5 site-selective samples. It is evi- 
dent that a remarkable oxygen isotope shift of T c as well 
as of cr sc is present. More importantly, the data points of 
the site-selective 16 O p 18 O ac ( 18 O p 16 O ac ) samples coincide 
with those of the 16 O pac ( 18 O pac ) samples. In order to put 
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FIG. 4 

versus the OIE exponent ao (-AT C /T C ) for 
Yo.6Pro.4Ba2Cu307_i samples. The error bars of the "triv- 
ial" 16 P ac point (zero isotope shift by definition) indicate 
the intrinsic uncertainty of the experimental procedures on 
the estimation of T c and \ a t- 
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TABLE I: Summary of the SOIE on T c and A~ 2 (0) for Y , 6Pro.4Ba 2 Cu 3 07_i extracted from the experimental data (see text 
for an explanation). The values of AT C /T C , ao, AA~ b 2 (0)/A~ i) 2 (0), and /3o are not corrected for the uncomplete oxygen-isotope 
exchange (see text). 



Sample 


T c 
(K) 


AT C 
(K) 


AT C 
T c 

(%) 


ao cr sc (0) 
(Ms- 1 ) 


Aa sc (0) 

(m s_1 ) 


(%) 


IbQ 


44.62(12) 






1.090(7) 






16q 18q 


44.66(10) 


0.04(16) 


0.09(36) 


-0.007(29) 1.070(9) 


-0.020(11) 


-1.8(1.0) 0.14(8) 


18q 


42.96(12) 


-1.66(17) 


-3.72(40) 


0.298(32) 1.022(9) 


-0.068(11) 


-6.2(1.0) 0.50(8) 


18q 16q 


43.17(11) 


-1.45(16) 


-3.25(36) 


0.260(29) 1.018(7) 


-0.072(10) 


-6.6(0.9) 0.53(7) 



are observed for the completely 18 substituted ( 18 O pac ) 
and the site selective 18 O p 16 O ac samples, which are the 
same within error bars. The other site-selective sample 
( 16 O p 18 O ac ) shows zero (within error bar) OIE on T c and 
a very small, but possibly non zero, OIE on A~ b 2 (0). This 
effect very likely arises from the small fraction of 18 O 
that substituted 16 O oxygen within the CuC>2 planes, as 
seen in the Raman data (see Fig. [I). However, within 
the error bars, the OIE on A~ 6 2 (0) can still be considered 
the same for the 16 O pac and 16 O p 18 O ac samples. 

In summary, we performed /iSR measurements of the 
site-selective oxygen isotope effect on the in-plane mag- 
netic field penetration depth in underdoped polycrys- 
talline Yo^Pro^E^CusOr-a samples. It was found that 
the whole or nearly the whole effect on the transition 
temperature T c and on the in-plane magnetic penetra- 
tion depth A^ 2 (about 100 % within error bars) comes 
from the oxygen within the superconducting CuOi planes 
and not from the apical and chain oxygen. Noting that 



the lattice parameters remain essentially unaffected by 
the isotope substitution^SiSi our results show the exis- 
tence of a strong coupling of the electronic subsystem to 
phonon modes involving movements of the oxygen atoms 
in the Cu02 plane, while suggesting that modes involving 
apical and chain oxygen are less strongly coupled to the 
electrons. This is in agreement with recent inelastic neu- 
tron scatterin g 10 ! 11 and photoemissioniSii^ studies, indi- 
cating a strong interaction between charge carriers and 
Cu-0 bond-stretching- type of phonons. 
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